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ABSTRACT 

Intermediate mass protostars, the bridge between the very common solar- hke protostars and the 
more massive, but rarer, O and B stars, can only be studied at high physical spatial resolutions in 
a handful of clouds. In this paper we present and analyze the continuum results from an observing 
campaign at the Submillimeter Array targeting two well-studied intermediate mass protostars in Orion, 
NGC 2071 and L1641 S3 MMS 1. The extended SMA (eSMA) probes structure at angular resolutions 
up to 0.2" , revealing protostellar disks on scales of ^200 AU. Continuum flux measurements on these 
scales indicate that a significant amount of mass, a few tens of Mq, are present. Envelope, stellar, and 
disk masses are derived using both compact, extended and eSMA configurations and compared against 
SED-fitting models. We hypothesize that fragmentation into three components occurred within NGC 
2071 at an early time, when the envelopes were less than 10% of their current masses, e.g. < 0.5 
Mq. No fragmentation occurred for L1641 S3 MMS 1. For NGC 2071 evidence is given that the bulk 
of the envelope material currently around each source was accreted after the initial fragmentation. 
In addition, about 30% of the total core mass is not yet associated to one of the three sources. A 
global accretion model is favored and a potential accretion history of NGC 2071 is presented. It is 
shown that the relatively low level of fragmentation in NGC 2071 was stified compared to the expected 
fragmentation from a Jean's argument. 

Subject headings: star formation, Orion, fragmentation, sub-millimeter interferometry, radio astron- 
omy 



1. INTRODUCTION 

Star formation research is a cornerstone of current- 
day galactic astronomy. A solid understanding of star 
formation allows us to analyze astronomical structures 
over a wide range of physical scales, from disk and 
(giant) planet formation around sunlike stars to the 
physical structure of the giant molecular clouds which 
form the building blocks of galaxies and which play 
an important role in galaxy evolution. Star forma- 
tion studies are necessary in order to define the initial 
conditions of most objects in our universe. Typically, 
galactic studies of star formation concentrate on either 

|tkempen@alma.cl| 



isokted low-mass (LM) (iLada fc Wilkins j 119841: iLada 
1987|: |Andre et all 119931: iShirlev et all 120 00: J0rgeiisen 



2OO4I fBottinelU et al.l 120041: lEvans et ai.lT2009i) or high 
mass (HM) star-forming regions, which are, with few 
exceptions, seen almost exclusively in clusters (e.g. 



Beuther et al.l2005l 120071 IRagan et al.ll200llSmith et al.l 
20091 : iKeto fc Zhang! \20m . Both ends of the mass 
spectrum provide unique perspectives on star forma- 
tion. Many isolated LM star-forming regions have the 
advantage of being nearby, allowing them to be ob- 
served at high physical resolution, and are far less af- 
fected by strong radiation fields created by neighbour- 
ing (proto)stars. These LM stars are often assumed t o 
form within individual collapsing envelopes (|Shulll977[ ). 
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Most field stars in our galaxy (see e.g iBressert et alj 
I2010D . however, form in clusters where interactions from 
nearby more massive stars, stronger radiation fields, 
and multiplicity due to fragmentation in t he parental 
cloud complica te the star f ormation proce s s ([Adams fc I 
[200lt Adams eTaD [200l iDuchene et all I2007D . Ah 
though physical models have existed for a decade (e.g. 
McKee fc Holhmanlll9 99l:lKlessenll200ll:lK rumhol z et al. 



2005; Bate & BonncU ~ |2005I: iVazguez-S cmadcn i et al. 
20091 and references therein), only a handful of re- 
cent studies have tried to observati o nally determine 
cluster properties ([Smith et al.l 120091 : [Longmore et al.l 
120111) . In broad terms, the evolution of isolated LM 
protostars is reas o nably well understood (iLadal 119871: 
Andre et al.' '1993"; 'Robitaillc ct al.l l2006t iCrapsi et all 
2008; van Kcmpcn ct al. 2009b). Even with the re- 
cent discovery of the VeLLO class, significant progress 
has been made to include these t ypes of protostars ir i 
the general evo l ution ary picture ([Dunham et al.l 120101 : 
iVorobvovl 120101 1201 1[ ). Recently an evolutionary pic- 
ture has been identified for very mass ive protostars 
([Fontani et all 120091: iKeto fc Zha ng"2010'). No system- 
atic effort, however, has been undertaken to observa- 
tionally confirm the theories of clustered star formation 
that haye been put forward by num erical modelling (e.g. 
IKlessenll200ll: IBate fc Bon]^l2005t [B ate 20 0l. 

Recent work from e.g. IKeto fc Zhang ([201 Of ) and 
iJohnston et al.l ([201 It ) show that there may be significant 
similarities in the formation mechanism of individual HM 
and LM protostars. As an example, for star formation 
at masses above M > 8 Mq there has been concern that 
radiation pressure might choke off the accumulation of 
mass from a surrounding envelope but recent results sug- 
gest that this is mitigated by the physical structure and 
geometry of the re gion, allowing O s tars to form throug h 
core collapse fe.g. lKrumh"oI3l2006t iKuiper et al.ll2010r ). 
The main observed differences in the star formation pro- 
cess all have an origin in the environment and the en- 
ergies involved in accretion and outflows. These can 
be enumerated as: (i) the clustering of LM protostars 
around HM protostars and thus the influence of the en- 
yironmental radiation field on LM protostellar formation 
([Krumholz et al.ll2010[) (ii) the fragmentation of natal en- 
velopes b efore nuclear fusion beg ins in the heaviest mem- 
bers fe.g. iBontemps et al.ll2010[ ) (iii) the strength of the 
internal radiation field produced by different stellar sur- 
face temperatures and its infiuence on the surroundings, 
and (iv) the feedback from powerful shocks created by 
outflow interactions with the parental cl oud, both induc- 
ing a nd dampening star formation (e.g lArce fc SargentI 

mm . 

Intermediate mass (IM) protostars (defined observa- 
tionally through their bolometric luminosity: Lboi> 50 
and < 2,000) have not been studied in detail, but may 
reveal crucial information on the differences between the 
modes of star formation. Although a few individual 
sources have been observed and analyzed in great de - 
tafl [e.g. NGC 7129 IRS 2, 'Fuente et alj ([2003 . 12001 . 
IRAS 20050+2720, Be hran e t al. (2008)], it is uncertain 
if these are either typical of their evolutionary stage/age 
or typical of clouds forming more massive stars than 
those in nearby star-forming regions, e.g. Taurus or 
Ophiuchus. Despite these limitations, IM protostars 
make excellent test-cases for star formation theories that 



aim to include th e full range of stellar masse s (e.; 
iFontani et alll2009l: Kama et al.ll2010l: iPalau et al.ll201C 
Being more luminous than LM protostars, these objects 
provide for larger warm zones within the enshrouding 
envelopes. These warm zones around forming stars were 
once thought to be nearly spherical in nature (the hot 
core), however, observations and models of key molec- 
ular lines now indicate that at least for LM protostars 
the structure of the inner en velope and outflow cavity 
walls plays an important role ("van Kemp en et al]l2009al : 
[Brudcrc r ct al. 2009; Visscr ct al. 2012). Recent obser- 
vations ([van Kempen et al] l2009aL l2010( ) show that this 
shell structure provides an accurate accounting of the 
heating and cooling balance and thus the evolution of 
low-mass protostars. Due to the large distance (~ 2 kpc) 
to most HM protostars, resolving the structure in such 
warm regions has proven difficult. IM protostars are nat- 
ural laboratories that can probe possible differences be- 
tween LM and HM protostars as they can be found at 
much smaller distances (0.5-1 kpc) than HM protostars 
and yet they produce significantly more UV photons than 
LM protostars. 

A second important characteristic of star formation 
theories can also be studied with IM protostars due to 
their proximity. Small mini-clusters of LM protostars 
cannot be distinguished from single HM protostars at 
large distances (> 2 kpc) with even the highest resolu- 
tions of the current generation of sub-mm interferome- 
ters. There are few observational papers investigating 
clustering and attempting to quantify the fragmentation 
of HM cores into multiple members e.g. iBeuther et al.l 
([200l : lBrogan et al.l ((20091 ) : lLongmoreetal.l([201 If ). IM 
protostars in Orion are near enough to allow facili- 
ties such as the SMA, the IRAM PdB interferometer, 
CARMA and in the future ALMA, to individually distin- 
guish cluster members and thus constitute an intriguing 
sample. 

This paper presents an analysis of the small-scale 
physical structure of two protostars in Orion, the near- 
est cloud producing protostars more luminous than 50 
Lq. The tw o observed sourc e s, NGC 2071 and L1641 
S3 MMS 1 (ISeth et al.l [200l IStoiimirovic"eraII i2008: 
iSkinner et al.l I20090 were selected specifically because 



(i) they are more massive than typica l low-mass pro- 
tostar s (e.g. the PROSAC survey in iJorgensen et al.l 
(2005", '2007) and pro tostellar surveys in Ophiuchus 
(Johnstone et al. 2000; ivan Kempen et al . 2009b)), (ii) 
they are located in the Orion cloud, where the ISRF 
is significantly larger than in other nearby r egions 
such as Taurus and Ophiuchus ([Jorgensen et al.| [2006). 
and (iii) they are still near enough (~ 450 pc) 
to allow the reasonable identification of circums tcUar 
disks (e.g.lM cCaughrean fc O 'dell' 19961: ISeth et al.1 120021: 
iMenten et al.. .20071 ) and protostellar envelopes. These 
two sources are not located in the most crowded region 
near Orion- KL, making it easier to separate the proto- 
stars from their surroundings. Further, these sources are 
thought to be representative of intermediate mass star 
formation in general and therefore excellent test cases to 
probe star formation in the mass range between the low 
and high-mass ends. In this paper, we present the con- 
tinuum observations and discuss the physical structure 
of these two sources. The format is as follows: § 2 dis- 
cusses the observations, the results are presented in § 3, 
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§ 4 analyses the data, and the structure of IM protostars 
is discussed in § 5. Finally, conclusions are presented in 
§6. 

2. OBSERVATIONS 

The two Orion sources investigated in this paper, 
NGC 2071 and L1641 S3 MMSl, were selected based 
on the proposed source list of the Herschel Guaran- 
teed Time key program WISH (Water in Star-forming 
region with Herschel, PI: E.F. van Dishoec k) subpro- 
grani on intermediate mass p r otostars (see e.g.lFich et alJ 
[2OT0t iJohnstone efaH \20m Ivan Dishoeck et all 1201 lE 
At a distance of d estimated around ^ 450 pc 
pohnstone et al.l 1200 ID , Orion is the only active star- 
forming region producing massive stars within 500 pc 
(jSadavov et al.ll2010l : iBuckle et ai]|2010[ ) and is thought 
to be the only cloud at such distances which is pro- 
ducing IM protostars. Note that the distances to indi- 
vidual regions of Orion might vary. The O MC might 
be as close as 414 pc (jMenten et all l2007f ). As such 
the error of this distance is likely as large as 40 pc . 
NGC 2071 is a muc h-studied area (jSneh fc Ballvl 119861: 
iTorrelles et al.lll998D and was very recently reported on 
bv iCarrasco-Gonzalez et all (|2012t ) using cm and 3 mm 
wavelength observations at very high resolution. The 
known properties of the two selected sources can be 
found in Table [TJ The two Orion sources were observed 
with the Submillimeter ArrajQ over three nights in Jan- 
uary/February 2010 using the compact and eSMA config- 
urations. Further observations of L1641 S3 MMSl were 
obtained in September 2010 using the extended configu- 
ration of the array. These SMA observations are comple- 
mented by continuum observations from 2MASS, Spitzer 
and the JCMT. 

2.1. SMA 

The Orion sources were observed in two campaigns 
on the SMA covering three configurations: compact, ex- 
tended, and eSMA. Table [2] describes the settings that 
were used during each of the observations, including 
beam sizes, bandwidths, correlator configurations as well 
as the bandpass, amplitude and gain calibrators used for 
the different dates. 

Initial observations were made in the compact configu- 
ration on 3'"^ and 4**^ January 2010. On February ll*'', 
observations were made in the eSMA configuration. In 
the eSMA configuratior0 the very extended configuration 
of the SMA is combined with the James Clerk Maxwell 
Telescope (JCMTJ3 and the Caltech SubmiUimeter Ob- 
servatory (CSO). Additionally, L1641 S3 MMSl was ob- 
served in a filler on September 28th 2010 using 6 dishes 
at 230.358 GHz in extended configuration. 

^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

2 The eSMA (extended SMA) is a collaboration of the SMA, 
JCMT and CSO, to join the three facilities into a single long base- 
lino sub-mm interferometer. F or more infor mation on the eSMA, 
see BottineUi et al. (2009) and lShinnaga et al. (2009) 

^ The JCMT is operated by The Joint Astronomy Center on be- 
half of the Science and Technology Facilities Council of the United 
Kingdom, the Netherlands Organization for Scientific Research, 
and the National Research Council of Canada. 



For the compact and extended configurations, instead 
of the normal 2 GHz correlator bandwidth, the enhanced 
double-bandwidth mode was employed. In this mode, 4 
GHz is obtained in both the lower and upper sideband for 
a total of 8 GHz bandwidth. The main molecular lines 
targeted were ^^CO 2-1, in the upper sideband (USB), 
and the two well-studied isotopologues ^^CO and C^^O, 
in the lower sideband (LSB). An analysis of these CO 
lines will be presented in a future paper. 

For NGC 2071, the quasars 0501-019 and 3C120 were 
used as gain calibrators and Uranus as a flux calibrator. 
For L1641 S3 MMS 1, the quasars 0607-085 and 0423-013 
were used as gain calibrators, and Titan as a flux calibra- 
tor. During the eSMA observations, the correlator setup 
was limited to the 2 GHz bandwidth mode due to the in- 
creased number of antennas. The limited observing time 
of 1.5 hours was nevertheless sufficient to cover a signif- 
icant part of the uv-plane. 3c273 was used as bandpass 
and Titan and Vesta as flux calibrators. Due to its close 
proximity, Vesta was also used as a gain calibrator. The 
correlator was fixed at 349 GHz, due to track-sharing 
with observations of Titan. 

The data reduction was performed using a combina- 
tion of the following software tools: the MIR package for 
IDL, MIRIAD, and CLASS in GILDA^. For all images, 
cleaning was done to 3a using the Clark method and a 
natural weighting was applied to obtain the best images 
for our goals. Although uniform weighting can provide 
better image fidelity, the price of the lower signal to noise 
and thus the detection of multiple components as well as 
their relative strengths, was considered not to be in line 
with the goals. 

2.2. Complementary data 

Several other astronomical data-sets are utilized in 
this paper as an aid to analyzing the source proper- 
ties. SCUBA 850 /xm continuum flux archive data 
(|Di Francesco et al.ll2008[ ) is available for NGC 2071. Un- 
fortunately, L1641 S3 MMS 1 was not included in the 
SCUBA archive due to its uncertain calibration. It was, 
however, observed by SCUBA and measurements were 
retrieved. iZavagno et all (|1997^ observed L1641 S3 MMS 
1 using the predecessor of SCUBA, the UKT14 common 
user bolometer instrument, and reported a calibrated 
flux of 5.14 Jy at 800 /xm. Comparing the UKT14 
results, corrected for the observed source size, and the 
uncalibrated SCUBA observations allow us to adopt a 
total flux of 5+/-1 Jy at 850^in for L1641 S3 MMSl 
and to calibrate the SCUBA map. Spitzer IRAC and 
MIPS 24 photometry were obtained for both sources from 
the Spitzer archive (Megeath, priv. comm). The data 
was taken within the scope o f the Orion Spitz er program 
[PI: Tom Megeath, see e.g. Gut ermuth et al.l (12009) and 
Megeath, T. et al. (in prep) for more information]. Near- 
IR photometric data points for the SED analysis were 
obtained from the 2MASS archive. The continuum mea- 
surements are provided in Tables |3] and HI 

3. CONTINUUM RESULTS 



GILDAS is a software package developed by 
IRAM to reduce and analyze astronomical data. 
|http://www.iram.fr/IRAMFR/GILDAS| 
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TABLE 1 
Source Properties 



Source 


R.A. 
hms 


Dec. 
dms 


^LSR 
km s~^ 


Lboi" 

L© 


Dist." 
pc. 


NGC 2071 
L1641 S3 MMl 


05 47 04.7 
05 39 56.1 


+00 21 44 
-07 30 28 


9.6 
5.3 


520 
70 


422 
465 



° Luminosities and distances arc adopted from the WISH list Ijvan Dishoeck et aLII20J^ and references therein. 







TABL 
Observational 


E 2 

PARAMETERS. 










Dates 


Config. 


Source 


Beam size 

("x") 


Bandwidth 
GHz 


Freq. 
GHz 


SB 


Band 


# Chan. 


Jan 3 2010 
Jan 4 2010 
Feb 11 2010 
Sept 28 2010 


Compact 
Compact 
eSMA 
Extended 


NGC 2071 
L1641 S3 MMSl 
NGC 2071 
L1641 S3 MMS 1 


3.4"x2.9" 
3.5"x3.0" 
0.37"x0.19' 
0.9"x0.87" 


4 GHz 
4 GHz 
2 GHz 
4 GHz 


230.538 
230.538 
349.415 
230.538 


USB 
USB 
USB 
USB 


sl3 
sl3 

sl7 
sl7 


128 
128 
128 
128 


Date 


Config 


Source 


Bandpass 


Flux Gal. 


Gain 


Calibrators 




Jan 3 2010 
Jan 4 2010 
Feb 11 2010 
Sep 28 2010 


Compact 
Compact 
eSMA 
Extended 


NGC 2071 
L1641 S3 MMS 1 
NGC 2071 
L1641 S3 MMS 1 


3c454.3 
3c273 
3c273 
3c454.3 


Uranus 
Titan 

Vesta, Titan 
Gallisto/Nept. 


0501-019, 3cl20 
0607-085,0423-013 

Vesta 
0607-085,0609-157 






™ w- 
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Fig. 1. — Comparison of continuum observations for both Orion sources across multiple resolutions. The first column shows the 850/im 
flux from SCUBA with a resolution of ~15". The second column shows the SMA-compact 1300/im observations, with a resolution of 
fs4.5",. The third column shows the very high resolution using the eSMA (NGC 2071 only) at a resolution of 0.3" at a frequency of 350 
GHz. Contour lines are shown at 0.1, 0.2... 0.9x the maximum flux at that resolution, as reported in Tablc[3] The synthesized beams are 
shown with black elipses, while the SCUBA beam is shown in the figure of NGC 2071. Extended configuration observations of L1641 S3 
MMS 1 can be found in Figure 4. 



For each of the Orion sources Figure [T] shows the sub- 
miUimeter continuum observations at different spatial 
resolutions, ranging from 15" to 0.2", using a combi- 
nation of the 850/im and 1.3 mm data from the SMA, 
eSMA and JCMT/SCUBA. The measured fluxes at each 
of these scales is tabulated in Tabled In the figure, im- 
ages are provided at three different length scales (120", 
20" and 2") and each image is normalized to the peak flux 
of the source in order to accurately compare the source 



structure. It is clear that the large-scale environments of 
L1641 and NGC 2071 look very similar at scales of '-SO" 
(Fig. [T]) and the peak fluxes of each core are remarkably 
close (5 versus 5.2 Jy), within the calibration error of 
SCUBA (20%). 

3.1. NGC 2071 

Multiplicity is observed with the NGC 2071 core when 
observed with higher spatial resolution using the SMA- 
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TABLE 3 

SUBMILLIMETER CONTINUUM MEASUREMENTS 



Source 



Core Offset 



Peak 
Jy/Bcam 



Total Flux 

Jy 



SCUBA 850 tim 



NGC 2071 
L 1641 S3 



5.8 +/- 0.8 
5^ +/- 1.0 



21.9 +/- 4.0 
30.0 +/- 6.0 



NGC 2071 



Compact 230 GHz 
A 4.2" 

B 1.4" 

C 12.2" 



L1641 S3 MMS 1 



1.4" 



0.45 +/- 0.1 5.6 +/- 1.1 

0.30 +/- 0.06 3.7 +/- 0.7 

0.14 +/- 0.04 1.7 +/- 0.4 

0.42 +/- 0.09 5.2 +/- 1.1 



Extended 230 GHz 






L1641 S3 MMS 1 1.6" 


0.26+/- 0.05 


1.8 +/- 0.4 


eSMA 349 GHz 






NGC 2071 A 4.2" 


0.13 +/- 0.04 


0.35 +/- 0.1 


B 1.4" 


0.12 +/- 0.04 


0.3 +/- 0.1 



1 The SCUBA flux L1641 S3 MMS 1 was estimated from the raw 
data. See text. Earlier papers (Snell &; Bally 1986; Trinid ad et al.l 



[20091 : [Uarrasco-Gonzalez et al.1120121'1 identify 2071-A with IRS-3, 
2071-B with IRS-1 and 2071-C with IRS 2. 




Fig. 2.— Zoomed image of the NGC 2071 protocluster at 230 
GHz continuum in the SMA compact configuration. Contour lines 
are in 3,6,9... a with a = 0.013 Jy/beam. Three cores are iden- 
tified, labelled A, B, and C and referred to in the text as 2071-A, 
2071-B, and 2071-C respectively. The phase center is indicated 
with a white plus sign. 

compact observations, as shown in the 2nd column of 
Figured! The source shows two fragmented peaks, 2071- 
A and 2071-B, with a third peak, 2071-C, to the north- 
east (see Fig. [5]), clearly detected (S/N>30) in the map 
which has a S-u rms noise of 4 mJy/beam. At the 
highest resolutions available (eSMA at '^840 /im) in- 
dividual disks in the NGC 2071 sources are discerned 
within the 2071-A and 2071-B cores. This eSMA map 
has a l-CT noise level of 7 mJy/beam. The error in the 
peak fluxes are dominated by the errors in flux cali- 
brations. These three peaks can be identified with older 
identifications (,Snell fc Bally, 1986: Torrcllcs et al. .. 19981: 



iCarrasco-Gonzalez et al.l 120121) : 2071-A is identified as 
IRS 3, 2 071-B as IRS 1 and 2071-C as IR S 2. The source 
shown in lCarrasco-Gonzalez et akl (l2012fl between A a nd 
B, first identified as VLA 1 bv iTrinidad et aD (|2009( ) is 
unresolved in the SMA map, although a non-circular 
elongation of sources A and B does point towards this 
position. In the eSMA map, a 3.4 a unresolved signal is 
seen at the position of VLA 1 . 

Figure [3] shows, for the cores 2G71-A and 2G71-B, the 
observed visibilities from the SMA and eSMA as a func- 
tion of projected baseline calculated from the flux peaks 
by using the UVAMP MIRIAD task. To derive these 
visibilities a gaussian model of the second source was 
subtracted in the image plane before calculating the vis- 
ibilities. The proflle for source A in the compact configu- 
ration shows a resolved envelope, up to 20 kA, combined 
with an unresolved central sour ce of 0.38 Jy (similar to 
Fig. 3 of iJgirgensen et all ()2005f )). Although source B is 
not as bright as source A, a similar fit can be reached 
in which an envelope is resolved up to ~30 kA, with an 
unresolved central component of 0.15 Jy. The eSMA ob- 
servations confirm that these unresolved components are 
disks, as they are resolved on scales of a few hundred 
kA corresponding to spatial scales of ~200 AU at the 
distance of Orion. The disks are both elongated, but 
suggest different major axis directions. Central unre- 
solved components within both 2071-A and 2071-B are 
constrained to < 0.1 Jy. 

If we extrapolate the 230 GHz fluxes to 345 GHz using 
a simple blackbody model with a /3 of 2, the total emis- 
sion picked up in the compact configuration SMA in A, 
B and C would be 19 Jy. With the observed flux of 21.6 
Jy in SCUBA, we retrieve about 90% of the total flux. 
Of course, this is highly dependent on the uncertainty of 
the absolute flux calibration. 

No disk-like component was detected down to a level of 
10 mJy/beam in the eSMA observations toward 2071-C. 
As this position was relatively far off the phase center 
(20 percent of the SMA antenna's primary beam) and in 
fact outside the primary beam of the JCMT and CSO, 
the presence of a disk c a nnot be ruled out. Note that 
ICarrasco-Gonzalez et al.l ()2012l ) resolve 2071-C as a close 
binary with a separation of about 500 AU. This would 
have been detectable with our resolution, but from their 
Fig 2 this would not be detectable with a 3 ct confidence 
with the sensitivity without the CSO and JCMT. We 
present a more detailed analysis in § 4.3. 

3.1.1. Spitzer maps 



Spitzer imaging from lSkinner et al.l (|2009[ ) can provide 
important constraints on the individual sources. Over- 
plotted SMA observations on the Spitzer-IRAC image 
are shown in Figure ID 

Comparing with the Spitzer data in Table [4] reveals 
that 2071-B and 2071-C are clearly detected in all bands, 
although 2071-B is dominated by emission at redder 
wavelengths. Source 2071-A, which is the brightest in 
the SMA images, is not clearly detected in any band, 
with only weak emission at 8 /im and 24 /xm. At the po- 
sition of 2071-C the source is dominated by the 3.6 /im 



emission. 
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Fig. 3. — Observed continuum visibi]ity amplitudes as a function of tlie projected baseline lengtli for 2071-A and 2071-B. Error bars are 
1 a statistical errors and tlie dotted histogram indicates the zero-expectation level at 1 a. Left column: The compact configuration data 
of sources 2071-A (up) and 2071-B (down). The data clearly indicate a classical profile as seen in e.g. Fig. 3 of j0rgensen c t al. (2003) 
where a resolved envelope transitions into an unresolved disk-like component. For 2071-A in the compact configuration this is indicated by 
a horizontal line. Richt column: The eSMA data for the two sources. Here it can be seen that the disks are resolved and the visibilities 
fall off until kA ~ 500, or about 80 AU at the distance of Orion. 




Fig. 4. — SMA continuum observations at 230 GHz at 3,6,9... a plotted in contours of the compact configuration in Fig. [T] and [2] over a 
Spitzer IRAC 3.6, 4.5 and 8 /xm false color image of NGC 2071. 
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TABLE 4 
Spitzer and 2MASS Photometry 



Source 


J 


H 


K 


IRAC 1 


IRAC 2 


IRAC 3 


IRAC 4 




Jy 


Jy 


Jy 


Jy 


Jy 


Jy 


Jy 


2071-A 


1.6C-4 


7.9e-4 


4.6e-3 


7.0e-2 


4.9e-2 


2.2e-l 


0.9 


2071-B 


3.0C-4 


2.2C-3 


2.1e-2 


4.4e-2 


3.7C-2 


1.7C-1 


0.7 


2071-C 


2.3C-3 


6.0C-3 


4.5e-2 


8.6e-2 


6.0e-2 


2.7C-1 


1.2 


L1641 S3 MMS 1 


- 


- 


- 


0.69 


0.82 


2.3 


17.2 
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Fig. 5. — The eSMA continuum observations at 349 GHz at a 
resolution of ~0.25" of the area around 2071-A and 2071-B, iden- 
tifying these cores, as well a s showing the locations of the peak 
fluxes at 3 mm as shown in ICarrasco-Gonzalez et aLl I I2012I) . In 
addition the Phase center is shown by a plus sign. Contours are in 
levels of 3-a-, with a equalling 7 nijy/beam. The synthesized beam 
is shown with a black ellipse. 

The recent paper by ICarrasco-Gonzalez et alj ()2012f ) 
identified for the fir st time the radio source VLA-1 
(jTrinidad et alJ [2Q09D at mm wavelengths, with a flux 
density of 9 mJy. At this position, a signal of 27 mJy, 
equalling a 3.85cr detection, can be seen, likely associated 
with VLA-1. Fig. [5] shows the eSMA image of the area 
around 2071-A, 2071-B and this potential new source . 
From the SED of Fig 2 in lCarrasco-Gonzalez erall (|2012[ ) 
a predicted flux of ~30 mJy can be derived at the eSMA 
wavelength. The detection at submillimeter would con- 
firm the theorized existence of another protostar in the 
NGC 2071 cluster. None of our other observations, in- 
cluding Spitzer, has the resolving power to distinguish 
any emission from VLA-1 from e i ther 2 071-A or B. X- 
ray emission from iSkinner et al.l ()2009[ ) has also been 
associated with this source. We therefore tentatively 
support the conclusions from ICarrasco-Gonzalez et aLl 
(J2012) that VLA-1 is a very young, even more embed- 
ded protostar. High-resolution, mid-IR observation are 
needed to distinguish any IR emission from the region to 
confirm this however. For the remainder of this paper, 
we will not discuss this source in any detail due to the 
limitations of the eSMA detection and observations at 
other wavelengths. Any other potential clumps similar 
to VLA-1 are not seen down to our detection limit. 

3.2. L1641 S3 MMS 1 

In contrast to NGC 2071, Figure [D reveals that L1641 
S3 MMS 1 is not fragmenting into multiple sources at 
small scales. Only a single continuum source is seen at 
a size scale of 4.5". The 1.3 mm continuum emission 
from the SMA compact configuration shows an almost 
perfectly circularly symmetric source. In extended con- 



FlG. 6. — Zoomed image of the L1641 protostar at 230 GHz con- 
tinuum in the SMA extended configuration. A plus sign indicates 
the phase center. Contour levels are at 10%, 20%, 30%, ... of the 
peak flux (0.26 Jy/bcam) 
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Fig. 7. — Observed continuum visibility amplitudes as a function 
of the projected baseline length for L1641 S3 MMS 1 in compact 
configurations. Error bars are 1 a statistical errors and the dot- 
ted histogram indicates the zero-expectation level. The extended 
configuration visibility was found to be unresolved at a level of 0.2 

Jy- 

figuration, Figure El the source looks unresolved , with a 
peak flux of 0.26 Jy/beam, roughly half of the peak seen 
in compact configuration of 0.45 Jy/beam. The central 
unresolved component is likely associated with 
disk with radius <300 AU while the more than 50% 
envelope is significantly filtered out by the extended con- 
figuration. It is also clear that the original coordinates 
of L1641 S3 MMS 1 are oflf by almost 2 arcseconds. The 
new updated position is 05'*39™56.1", -07''30"28". The 
original coordinates were based on the IRAS catalogue 
with a 2 arcminute beam. Radio continuu m observations 
strong enough to be used for astrometry (jMorgan et al.l 
,1990.) also identify the radio continuum source within < 
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1" of our position. 

The visibihties shown in Figure [7] confirm that the en- 
velope is resolved and that the structure variations are 
extend inwards to 60 kA, a significant difference with 
NGC 2071. No unresolved component dominates the 
emission between 20 and 60 kA. Any unresolved com- 
ponent that could correspond to a central disk is limited 
to <0.1 Jy and <400 AU in size. The envelope structure, 
as seen by resolved out emission at shorter baselines, ex- 
tends all the way inwards to 700 AU. 

The Spitzer image for L1641 S3 MMS 1, Figure El 
shows that the embedded source is associated with the 
bright greenest emission (4.5 /Ltm). The 4.5 micron IRAC 
band contains spectral features which emit stro ngly in 
shocked molecular gas (|De Buizer fc Vaccall2010[ ). so ex- 
cess emission in this band is ofte i i attributed to outfiows 
from YSOs (jCvganowski et al.l 120081 : iChambcr s et al.l 
|2009| ). Other potential members in the L1641 protoclus- 
ter can be seen, but none emit at 1.3 mm wavelengths. 
Megeath et al. (in prep) classifies these other sources as 
normal stars or TT auri stars. The flow was earlier also 
seen at H2 emission (jStanke et al.ll2000l ). 

4. ANALYSIS 

4.1. Mass derivations 

If one assumes that the cold dust has the emis- 
sion properties proposed b y column 5 of the table 
lOssenkopf fc Henninal (|1994[ ). also known as 0H5 dust, 
one can calculate the observed mass at each physical scale 
shown in Figure [TJ The mass can be calculated directly 
from the continuum emissio n using the following formula 
(see e.g. iShirlev et ani2000[) : 



TABLE 5 

Mass estimates 



Md = S,DyB,{Tdust)n,, 



(1) 



where S^, is the integrated flux density, Bi,{Tdust) is the 
Planck function at the dust temperature Tdust, and Ki, is 
the opacity per gram of gas and dust. It is assumed 
that at both observed wavelengths (1.3 mm and 850 
fim) the dust is optically thin. For 0H5 dust, k^, 
is 0.02 cm^ gm-i at 850 fim and 0.009 cm^ gm"! at 1.3 
mm. Without additional measurements, however, the 
dust temperature cannot be independently derived. At 
large scales, the mass is dominated by the outer cold 
du st and the temperat ure is typically assumed to be 20 
K (jShirlev et al.|[2000l ) and we adopt that value in this 
paper. Material heated to higher temperatures, e.g. by 
outflow shocks, will also emit at these wavelengths but 
should not account for a large fraction of the mass. Thus, 
we assume that the bulk of the material is dominated by 
the cold gas on all spatial scales. Table [5] provides the 
derived masses as a function of spatial scale. 

The derived total mass is likely to have a significant 
level of uncertainty due to uncertainties in the observed 
fluxes, the gas to dust ratio, and the underlying range of 
possible dust properties. We estimate the uncertainty in 
the total mas s at about 40% from t hese effects and refer 
the reader to IShirlev et all pOOOl ). IShirlev et all pOOl 

^ Note that th i s ma y be a source of error. E.g. 

[Andrews &: WiUiamsl 1)20051 ') found that about 25% of the emission 
of dust may be optically thick. It follows that even at 1.3 mm, there 
is a sizeable contribution of the optically thick emission. However, 
for simplicity we continue with the assumption that the emission 
is optically thin. 



Source 



15" 4"" 0.9"° 
Mass [Me] 



0.25 ' 



NGC 2071 



Core" 21.7 

A 

B 

C 



12.3 

8.2 
3.7 



0.35 
0.29 
<0.2 



L1641 S3 MMSl 



20.9 11.46 3.9 



" Observations taken at 230 GHz with dust opacity k=0.009 cm 

gm"^ 

^ Mass measurement derived for the single protocluster core. 



and IDi Francesco et al.l (j2008[) for a more in-depth dis- 
cussion. The uncertainty in the dust temperature, how- 
ever, can introduce a significant uncertainty in the de- 
rived mass. Using a lower dust temperature of 10 K 
increases the mass by a factor of about three, while as- 
suming a higher temperature of 30 K yields about 2.5 
times less mass. 

On large scales the masses are dominated by the cold 
dusty envelopes. Only small fractions (<0.5 Mq) re- 
main visible at the highest resolutions observed with the 
eSMA. Even if a large amount, e.g. 75%, of the disk 
material is resolved out by the lack of shorter baselines, 
disk masses would still be below 2 M0. A scenario in 
which the disks are largely resolved out is unlikely, since 
at the distance of Orion, ^^450 pc, the resolved eSMA 
disks are still ^200 AU in radius. Within the extended 
observations of L1641 S3 MMSl, which probes down to 
scales of 400 AU, the mass ('^3.8 M0) is still dominated 
by the cold outer envelope, extending down to a few hun- 
dred AU. An unresolved disk( < 400 AU) would have a 
mass <0.45 Mq, assuming the unresolved inner part has 
brightness ^0.2 Jy as seen from the longest baselines. In 
conclusion, all sources, independent of possible IR detec- 
tions with Spitzer are clearly deeply embedded sources 
equivalent to the Stage 0. 

4.2. Spectral Energy Distributions 

A powerful probe of protostellar evolution is the 
Spectral Energy Distribu tion (SEP). St u dies, such 
as demonstrated in e.g. iRobitaille et al.l ()2007[ ) and 
lEvans et all (|2009[ ). have used databases of 2MASS and 
Spitzer sources to efficiently characterize large number 
of protostars using their SEDs. Although som e stud- 
ies, such as iCrapsi et al. (2008), lEnoch et all UdOS), 
Ivan Kempen et all (j2009b ). and Maxwell & Johnstone 
(in prep) , have shown that individual SEDs may not pro- 
duce unique classification and physical structure param- 
eters, they nevertheless do provide a reasonable guess 
at the evolutionary stage and structure of the central 
protostar. Using the Spitzer data, complemented by 
2MASS data where available, SEDs in the range of 1 
to 100 /im were constructed for all four detected sources 
(2071-A,B and C and L1641 S3 MMS 1) and in turn 
fitted with a model from the SED grid of models pre- 
sented by Robitaille et al.l ([2006) using the online SED 
fitteiQ (jRobitaille et al.ll2007[) . The submiUimeter obser- 
vations were not included due to either their uncertain 

^ The SED fitter can be found at 

|http://caravan.astro.wisc.edu/protostars/| 
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Fig. 8. — SMA continuum observations at 230 GHz at 3,6,9... a of the compact configuration Fig. [T] plotted in contours with a = 2.3 
mjy/beam in over a Spitzer IRAC 3.6, 4.5 and 8 ^m false color image of L1641 S3 MMSl 
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Fig. 9.— SED fits for {Top Left) L1641 S3 MMS 1, {Top Right) 2071-A, {Bottom Left) 2071-B, and {Bottom Right) 2071-C. The black 
line shows the best-fitting model, the grey lines show the next 10 best fits. Solid circles show confirmed fiuxes, while triangles show either 
lower (normal triangle) or upper limits (upside-down triangle). Bars show the uncertainty in the fiux determination. The dashed line is 
the stellar SED if it would be visible unobstructed by circumstellar material. 
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calibration (SCUBA) or the lack of short spacings and 
the filtering out of emission at larger scales (SMA) . 

Due to the possi ble i naccu racy of individual SED fits 
shown by CraDsi e t alJ (120081) . as well as the duplicity of 
many fits ( Robitail le'et al.ll2006ll2007D . we average the 10 
best SED results and study only the emission shortwards 
of 100 /im(models which grossly under- or overestimate 
the long wavelength fluxes are dismissed before the av- 
eraging). Figure 14.21 shows the fitted SEDs for L1641 
S3 MMSl, 2071-A, 2071-B and 2071-C. Table H in turn 
presents the resulting envelope masses, bolometric lu- 
minosities, stellar masses and the ratio of envelope over 
stellar mass of the SED, derived by averaging the 10 best 
fits. 

Table |6] reveals that the combined luminosities of the 
three NGC 2071 sources do not account for the pre- 
viously derived bolometric luminosity of 520, assumed 
for the total core. Most of the high core luminosity on 
larger scales is due to the inclusion of IRAS 60 and 100 
/imphotometry, which includes a significant contribution 
from the ISM and cloud material surrounding the core. 
Even the removal of these points indicates that the bolo- 
metric luminosity is dominated by large-scale cold cloud 
material. Higher resolution imaging using the PACS and 
SPIRE instruments on Herschel, as proposed by HOPS 
(PI: T. Megeath) will re-observe these regions and deter- 
mine more accurate core luminosities. 

The envelope masses found via SED fitting and direct 
measurement of dust continuum emission range from a 
few solar masses for 2071-C to almost 15 Motor 2071- 
B. This is significantly higher than the masses com- 
monly found in ne arby clouds, which are typically a 
few tenths t o -2 M^ HShirlev et al.ll2000H Johnstone et al.l 
120001 I200TI ). However, the masses derived from the SED 
are within 20 % of the masses derived from the compact 
SMA emission at 230 GHz. 

Although the Robitaille SED fitter can be used to de- 
termine limits on the disk emission, the high Ay and in- 
creasing optical depth of the dust makes this very uncer- 
tain and untrustworthy for these observations. We do, 
however, derive estimates on the stellar masses. From 
the SED fitting, the L1641 S3 MMS 1 internal source is 
significantly more massive (3.5 M©) than any of those 
in the NGC 2071 minicluster (< 1 Mq) and the ratio 
of envelope over stellar mass is much lower (this is also 
true if the SCUBA measurements are used to determine 
the envelope mass). This may indicate that L1641 S3 
MMS 1 has had a different accretion history than the 
NGC 2071 cluster members. Either it is older, or the 
accretion rate was significantly higher. Changes in the 
accretion history produce a significant influence on the 
radiation fleld produced by the star at UV wavelengths. 
We will be investigating this effect using molecular line 
spectroscopy in a forthcoming paper. 

5. DISCUSSION ON FRAGMENTATION 

Although a few individual IM sources have been 
observed an d analy z ed in detail (e.g. NGC 7129 
IRS 2. iFuente et all (l2005l 120071) . IRAS 20050-1-2720, 
IBeltran et al.l ()2008[ )). it is uncertain if these are (i) typi- 
cal of their evolutionary stage or age, (ii) typical of clouds 
forming more massive stars than those in nearby star- 
forming regions, e.g. Taurus or Ophiuchus (iii) typical 
for fragmentation at small (~ 500 AU) scales. Ideally, a 



much larger sample of IM sources needs to be considered, 
and these objects must be investigated at higher spatial 
resolution. 

The interferometric observations of NGC 2071 and 
L1641 S3 MMS 1 allow for an investigation into their 
fragmentation histories. Although similar on large scales, 
these two IM sources are quite different on small scales, 
with one separating into multiple components and the 
other remaining a single protostar. In this section we 
will attempt to unravel the fragmentation history of two 
cores massive enough to form an intermediate mass pro- 
tostar and determine which physical parameters make 
this possible. For an extensive in-depth discussion of 
fragmentation, we refer the reader to the s eries of papers 
by Alan Boss (e.g. 'Bosf^'1993', "2005, "2009^ and Matthew 
Bate (e.g.^atc 1998; Bate & Bonnch 2005.: .Bate. 2009. ). 

5.1. Conceptual Model 

To make the problem of molecular cloud fragmentation 
tractable, we start by assuming the intrinsically contin- 
uous fragmentation process of an original core can be 
separated into N discrete time steps: t=to, ti, t2,... end- 
ing with what is observed today at t=tnow Assuming the 
gas is of uniform density, isothermal and thermally sup- 
ported against gravitational contraction, the mass and 
separation of gravitationally induced fragments is given 
by the Jeans equations, written in terms of density and 
temperature. 



M, 



( T 



VlOK 



(10?^)"^^ 



(2) 



and 



A , = 2 X 10"* 



T 
TOK 



(toI^)"'^^ (^) 



It is important to note that both Mj and Aj are min- 
imum requirements for gravitationally induced fragmen- 
tation and that for spherical condensations without inho- 
mogeneities the fastest growing collapse mode consists of 
the entire core. For non-spherical condensations, or cores 
with significant inhomogeneities the opportunity for frag- 
mentati on at these Jeans scales is much enhanced (see for 
example iPon et al.|[20ll[) . 

The density, n, and temperature, T, of the gas will 
change with time due to the collapse itself, as well as 
the additional infall of large-scale material. The average 
density and temperature at each time-step can be de- 
noted To and ng, Ti and rii, ... Depending on the local 
physical conditions at a given time the cloud may break 
up into ever smaller fragments. We denote the observed 
fragments as A,B,C, .... and the precursor as the sum 
of final fragments (i.e. for three observed fragments A, 
B, and C, the initial core would be designated ABC). 

From the given parameters, one can derive two basic 
scenarios that a cloud can follow during fragmentation 
: (1) Direct fragmentation: A cloud directly fragments 
into the resulting structure in a single time-step. In 
this scenario the conditions for fragmentation are reached 
throughout the cloud and the resulting cores are all co- 
eval. There is only a single event ti. Notation: ABC 
=> A+B+C; (2) Hierarchical fragmentation. In this sce- 
nario it is possible for a fragment created in time-step tx 
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TABLE 6 

Results from Robitaille SED fit" 



Source 


Lboi 


Env. 


Mass*" 


Stellar Mass 


(Menv/Mstar)'' 




Lo 


Mq 




M© 




2071-A 


<27 


8.2 




0.9 


30 


2071-B 


10 


14.2 




0.5 


20 


2071-C 


3.4 


~3.5 




0.5 


6.8 


L1641 S3 MMS 1 


250 


>9.6 




>3.5 


2.7 



" Results are averaged over the 10 best-fitting models. The spread in these values is about 50% from the given value. Ay-values range from 
5-35 and are the main cause for the large uncertainty. 
^ Lower limit estimates. 



to fragment again into smaller fragments at a later time 
steps, which we label t^^i, ta;2,.... This can be done as 
long as the requirements as defined by the Jean's equa- 
tions are met. Notation: ABC =^ AB+ C or AC+B or 
A+BC^ A+B+C. 

A fragmentation event is characterized by two num- 
bers, Aj and Mj. The number of possible fragments 
created can be determined by dividing the mass of the 
cloud by Mj, however, as noted above the actual de- 
gree of fragmentation is determined by inhomogeneities 
in density and geometry (assuming for the moment a 
cloud held up only by thermal support and thus neglect- 
ing additional support mechanisms such as rotation or 
magnetic fields). The final mass of a fragment is likely 
to be much higher than the instantaneous Jeans mass, 
both due to this under-fragmentation effect and the ad- 
dition of further material from the accreting core. The 
fragmentation scale, however, is expected to provide a 
more durable result - close fragments must form at high 
densities and low temperatures. 

5.2. Applying the conceptual model to NGC2071 

Returning to the case of NGC 2071, we observe three 
present-day key components: A, B and C (see Fig 5). Al- 
though a tentative detection of a fourth source, VLA-1 
(or D), is discussed above, the lack of emission at infra- 
red wavelengths uniquely tied to this position means we 
cannot rule out this being a shock position. The applica- 
tion of the conceptual model for three or four fragments 
is identical, including the caveats discussed later. We 
assume these formed from the same initial core, ABC, 
at t=to. This means there can only be either one or 
two fragmentation steps: either i) ABC directly frag- 
mented into A, B and C, or ii) ABC first fragmented in 
two (AB-I-C, AC-I-B or A-I-BC), and one of these subse- 
quently fragmented. We further assume that the sepa- 
ration of the cores, Axyi has not changed significantly 
between the fragmentation event and the time of obser- 
vation: fobs- 

We first note that the fragments appear to have a hi- 
erarchical structure, with fragments A and B, near the 
core center, with a small separation Aab ~ 2500 AU and 
fragment C, closer to the core edge, with a separation 
^{AB)c ^ 5000 AU from the AB pairQ Considering the 
Jeans equations above, the simplest way to account for 

"^ The formation and fragmentation of a mini-cluster is formed in 
three dimensions and not in two. As such the observed distances, 
and the corresponding Jeans lengths, A, are projected distances 
and not physical distances. The difference in projected and phys- 
ical distance between AB, BC and AC can of course be different. 
If one assumes there is a significant depth in the core , distances 
can be assumed to be larger by a factor l/sin{i) with i the angle 



the separation variation is to assume that the density 
was higher by a factor of four (or the temperature was 
lower by a factor of four) when A and B fragmented com- 
pared with C. Similarly, the best handle on the density at 
fragmentation is supplied by the separation of the frag- 
mentcl, as noted in the previous section, and thus we 
suggest that the fragments formed when the gas density 
was at least 3 x 10^ cm~'^ (assuming an isothermal gas at 
T = 20 K). For these conditions, the instantaneous Jeans 
mass of each fragment would only be ^ 0.5 Mq, suggest- 
ing that additional accretion of material onto each frag- 
ment was required. To reproduce the higher masses 
of the central fragments A and B, compared with frag- 
ment C, utilizing only direct fragmentation requires ei- 
ther lower initial densities or higher temperatures. Con- 
versely, the smaller separation between the fragments, 
and their proximity near the center of the core, suggest 
higher density and lower temperature conditions. The 
simplest explanation for the mass discrepancy is that the 
original core did not fully fragment into M/Mj pieces 
(where Mj r^ 0.5 Mq), and thus the masses of A and B 
do not reflect the instantaneous Jeans mass but rather 
reflect the bulk mass of the inner core after further ac- 
cretion. Note it is possible D is one such fragment, but 
that due to its location it was deprived of incoming mass 
by A and B. Under such a scenario, the original core 
would have had only mild inhomogeneities, resulting in 
only the few fragments observed. One of these fragments, 
relatively far from the core center accreted little material 
beyond its initial fragment whereas the two central frag- 
ments continued to accrete the bulk of the core material. 
We thus propose the following cartoon history for NGC 
2071, as shown in Fig. [TOl The original supercore, ABC, 
collapsed isothermally until reaching a critical density 
ni > 3 X 10^ cm~^ in the central region (~ 5,000 AU) 
at ti where it fragmented into the cores AB and C with 
the observed separation. At t2 the mean density had 
increased by a factor of four and AB fragmented into A 
and B. After t2. A, B and C formed protostars, accreting 
further mass from the natal supercore. The heating of 
these central stars inhibited further fragmentation due to 

between two cores along the line of sight. On average this will be 
30% and at this level of uncertainty these results still hold. 

* It is of course possible that there is a dynamical component 
that changes the location of either A, B or C after formation. From 
unpublished data we observe a small difference of <0.15 km s~^ 
between C and AB, while there is no detectable difference <0.1 km 
s~^ between A and B in line observations of C^*0 (van Kempen et 
al. in prep). This limits the movement of C as compared to AB to 
<3,000 AU, assuming a lifetime of 0.08 MYr (50% of the average 
Class lifetime (Evans et al. 2009)). It is therefor not possible to 
rule out that Core C was ejected by A and B, although it is unlikely 
given the calibration errors in the velocity determination. 
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Fig. 10. — Cartoon representation of the fragmentation of NGC 2071. Masses are in units of Mg 



the increased gas temperature, since the Jeans' lengths 
of remnant density irregularities increase to scales larger 
than the core. During all this the supercore may well 
have continued to accrete from the global mass reservoir 
seen at larger scales. 

Whether the core fragmented hierarchically in time is 
debatable. It is likely that the core initially fragmented 
into AB and C with a later fragmentation of AB into 
A and B when the central density had increased signif- 
icantly. However, observed pre-stellar cores are known 
to have increasing dens ity toward their cen ters (and de- 
creasing temperatures (jBergin et al.ll2006() ) and thus a 
single fragmentation scenario cannot be ruled out. 

In the above analysis we have inherently assumed that 
the thermal pressure is dominant. However, based on 
single dish measurements and our SMA molecular line 
data (van Kempen et al. in prep.), we know that cores 
are highly supersonic {i 5). Clearly the cores are not 
thermally supported. Including such a significant non- 
thermal motion as an equivalent pressure in eqn. (2) and 
(3) of Jeans analysis has been shown to be more consis- 
tent with the observed core masses and separations (e.g. 
iPillaiet al.|[20lll) . From our SMA C^^^O data, the line 
FWHM is 2.5 km/s, corresponding to a velocity disper- 
sion of 1.1 km /s. For a density o f 10^ cm~^ following eqn. 
(3) and (4) of ' Pillai et all (|2011[ ) , we then derive a Jeans 
mass of 100 Mq, much higher than the total mass of the 
clump as derived from archival SCUBA data. The corre- 
sponding Jeans length also would require the cores to be 
separated by >20,000 AU inconsistent with the observed 
projected separations. However increasing the density to 
few times 10^ cm~^ would give Jeans mass and length 
consistent with the observed values. Such initial core 
densities are extreme (and unlikely) and therefore, we 
conclude that turbulent jeans fragmentation to present 
values would not justify our observations. 

Perhaps the most interesting aspect of this analysis is 
the suggestion that the NGC2071 core underwent much 
less fragmentation than could have been possible given a 
simple Jeans stability argument. The present core den- 
sity implies that the number of possible fragments is in 
the hundreds and even assuming somewhat lower initial 



density conditions, the number of fragments could have 
numbered more than ten. This suggests that there may 
be additional support mechanisms working against grav- 
itational fragmentation. Examples of support mecha- 
nisms include: rotation, turbulence, and magnetic fields. 
We will discuss these in the next section, as these support 
mechanisms are also the reason L1641 S3 MMS 1 did not 
fragment at all. Even with the inclusion of VLA-1 as a 
primordial fragment, these support mechanisms would 
still have prevented fragmentation at similar level. 

5.3. Why did LI64I not fragment? 

Given the large-scale similarities between NGC 2071 
and L1641 S3 MMS 1 it is curious that the later did 
not fragment. However, as noted at the end of the pre- 
vious section, even NGC2071 seems to have undergone 
less fragmentation than might be expected from a Jeans 
gravitational instability argument. From this perspec- 
tive it would seem that NGC 2071 and L1641 S3 MMSl 
are more similar than they are different and that in both 
cases support mechanisms against collapse are required. 
One additional parameter can also be raised for L1641 
S3 MMS 1 and that is the lack of inhomogeneity. 

We note that for low mass pre-stellar cores very lit- 
tle initial inhomogeneity has been observed in the den- 
sity stru cture beyond that required to support against 
gravity (jSchnee et al.l 120101 ). suggesting that these ob- 
jects will at least begin their collapse to protostars mono- 
lithically. The evolution of isothermal structures with 
constant density depends explicitly on the detailed ge- 
ometry as global, large-scale collapse, competes against 
gravitational fragmentation (See Pon et_al..2011, for a 
discussion) In general the more spherical the initial core, 
the less fragmentation expected. 

5.3.1. Support Mechanisms against fragmentation 

Additional support mechanisms against fragmentation 
which should be considered include rotation, turbulence, 
and magnetic fields . For example, 3D hydrodynamical 
simulations bv lBatd ()1998[ ) show that sufficient rotational 
energy in a cloud or core prevents fragmentation as the 
growth of the non-axisymmetric perturbations is stifled. 
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Gravitational torques actually remove the angular mo- 
mentum. Close binarie s can still be formed. Indeed 
iCarrasco-Gonzalez et al.l (|2012[ ) reveals that 2071-C is 

likely a close binary (~500 AU). 

Turbulence has a similar stifling effect (Bate 2009) by 
preventing or delaying the growth of non-axisymmetric 
perturbations. Turbulence can be introduced by e.g. am- 
bipolar diffusion of external radiation sources or the for- 
mation of the first protostar and associated outflow. The 
inclusion of magnetic fields highly complicates the forma- 
tion and fragmentation of protostellar cores in molecu- 
lar clouds. Similar to turbulence and rotation, magnetic 
fields are able dampen the fragmentation through the 
process of magnetic bra.king as well as magnetic p ressure 
(jHennebelle fc Tevssieii [20081: iPrice fc Batell2007[ ). But 
magn etic fields can also enhance fragmentation (jBossI 
12009^ . Magnetic tension prevents accretion and thus the 
fast growth of a central density singularity. This tension 
drives small density variations toward the Jeans' mass. 
A key parameter that determines the scale of the mag- 
netic tension and raking seems to be the initial shape 
of the molecular cloud with respect to the direction of 
the magnetic field. For instance, an oblate cloud will 
fragment much more than a prolate cloud. The study 
of fragmentation of detailed MHD codes is a very active 
field with many parameters that need to be included. A 
main conclusion is that the effect of magnetic fields is 
highly dependent on the initial conditions, in particular 
the small-scale structure and the inclusion of all physical 
characteristics, e.g. ambipolar diffusion. 

6. CONCLUSIONS 

In this paper we have presented millimeter and submil- 
limeter interferometry observations of two intermediate 
mass (IM) protostellar cores in Orion, L1641 S3 MMS 1 
and NGC 2071. The following conclusions can be drawn 
from the data: 

• Continuum observations from the JCMT/SCUBA 
and the SMA reveal that at large scales these two 
sources are similar in mass ('--^25 M0) but that at 
smaller spatial scales NGC 2071 has fragmented 
into three, or potentially four, low mass (LM) pro- 
tostars, which is not observed in L1641 S3 MMSl. 
Central stellar masses derived from SED fitting re- 
veal that NGC 2071 contains low-mass sources, 
with stellar masses < 1 Mq. The central star in 
L1641 S3 MMS 1 is > 3 M©. Ahhough it is possible 
the protostars in NGC 2071 will still form A-stars, 
only L1641 S3 MMSl will likely form a B7-B9 star. 
All three protostars in NGC 2071 remain deeply 
embedded, based on the the ratio of envelope over 
stellar mass, and thus Class or Stage sources a s 
defined by the classification of lEvans et al.l ()2009[ ). 
L1641 S3 MMSl is a much more massive protostar 
(stellar mass ^ 3 Mq) and has not fragmented. It 
also classified as Stage 0, despite being bright in 
the infrared (Spitzer). 

• Disks around the three sources in NGC 2071 are 



constrained to masses smaller than 0.35 Mq, about 
2% of the total mass of the envelopes and 10% 
of the derived stellar masses. The possible disk 
around L1641 S3 MMSl cannot be fully character- 
ized due to a lack of very long baselines, but is also 
constrained to < 0.45 M©, < 4 % of the total mass. 

• Calculation of Jean's masses and lengths and sub- 
sequent comparison to the physical distribution 
shows that much less accretion must have taken 
place on the NGC 2071 cluster as well as on L1641 
S3 MMS 1. In the latter case the stifling pre- 
vented fragmentation altogether. The stifling likely 
was caused by a combination of support mecha- 
nisms against fragmentation, such as rotation, tur- 
bulence and magnetic breaking. Fragmentation at 
very small scales (< 100 AU) cannot be ruled out. 

• The bulk of accretion onto the three NGC 2071 
cores took place after the fragmentation. These re- 
sults favor a model of competitive global accretion 
from a large-scale molecular cloud. Fragmentation 
might thus have taken place when each core was 
^0.5 Mq, about 10% of their current masses. 

Future studies of these sources should focus on observa- 
tions using even longer baselines and thus higher spatial 
resolution. Studies with e.g. the Atacama Large Mil- 
limeter/Submillimeter Array (ALMA) can probe down 
to resolution of 0.2-0.1 ". These will be able accurately 
characterize the protostellar disks surrounding the stars 
at high precision. 
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